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Defects in LiFeAs are studied by scanning tunneling microscopy (STM) and spectroscopy (STS). 
Topographic images of the five predominant defects allow the identification of their position within 
the lattice. The most commonly observed defect is associated with an Fe site and does not break 
the local lattice symmetry, exhibiting a bound state near the edge of the smaller gap in this multi- 
gap superconductor. Three other common defects, including one also on an Fe site, are observed 
to break local lattice symmetry and are pair-breaking indicated by clear in-gap bound states, in 
addition to states near the smaller gap edge. STS maps reveal complex, extended real-space bound 
state patterns, including one with a chiral distribution of the local density of states (LDOS). The 
multiple bound state resonances observed within the gaps and at the inner gap edge are consistent 
with theoretical predictions for s^ gap symmetry proposed for LiFeAs and other iron pnictides. 
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Impurity physics plays a key role in superconducting 
systems, beginning with the remarkable feature that non- 
magnetic defects do not strongly impact superconductiv- 
ity in conventional s-wave materials^. In contrast to such 
single sign s-wave superconductors, where only magnetic 
defects cause pair-breaking and in-gap states^, both po- 
tential and magnetic defects can induce in-gap states in 
d-wav^^H^ and multi-band sign reversal s-wave supercon- 
ductors (s^^PM^. Not surprisingly then, a superconduc- 
tor's sensitivity to defects, plus the energetic and spatial 
characterization of bound states localized at defect sites 
have provided clues to the pairing symmetry of novel su- 
perconductors, and STM has proven an invaluable tool 
for such studies^ ^ lo ii Since a sign change of the or- 
der parameter gives rise to sensitivity to defects, it has 
been suggested that the study of these impurity bound 
states in the iron arsenides could help close the on-going 
discussion regarding the gap structure (s++ or s^)'' ^. 

The fact that most high temperature superconductors, 
cuprates as well as pnictides, require chemical substi- 
tutions to tune them into their superconducting states 
adds further importance to understanding the effect of 
defects in these systems^^. In the cuprates, this tun- 
ing is often achieved through cation substitution on 
sites away from the Cu02 planes, doping them with 
holes or electrons, while largely avoiding strong scat- 
tering. Direct substitution onto the Cu02 planes is 
typically pair-breaking, sometimes strongly, sometimes 
weakly^. For many iron-based superconductors chemi- 
cal substitution that suppresses extended magnetic order 
is an essential ingredient in achieving high temperature 
superconductivit}/^^. In BaFe2As2 it has been shown that 
certain elements such as Co and Ni substituted into the 
Fe-layer induce superconductivit}'^^ while other sub- 
stituents (Mn) cause strong pair-breaking^A Thus, the 
arsenides lack the easy distinction of off-plane substitu- 



tion to promote superconductivity versus on-plane de- 
fects that are pair-breaking. This makes it particularly 
important to assess individual impurities and their influ- 
ence on the surrounding electronic states. 

STM study of the pnictides has proven difficult due to 
surface specific effects arising from a lack of natural cleav- 
ing planes, or from structural or electronic reconstruction 
caused by a polar catastrophe^. The issue has been re- 
cently solved with the discovery of suitable surfaces in 
cleaved LiFeAs crystal^^, and molecular beam epitaxy 
grown films of FeSd^ and KFe2Se^. All three also 
possess the advantage of being superconducting without 
chemical substitutions. In LiFeAs, STM has been used 
to m easur e the superconducting gaps of clean defect free 
area d^^ * ^^ , to study vortice^^, and to determine band 
structure and gap symmetry through quasiparticle in- 
terference induced by defect scattering^^ A detailed 
investigation of the impurities themselves and their local- 
ized electronic effects has yet to be reported, though in- 
gap states have been observed for iron adatoms in FeS^^ 
and possibly iron vacancies in KFe2Se^. 

In this letter we characterize defects arising from crys- 
tal growth in nominally stoichiometric LiFeAs. We iden- 
tify their position in the crystal lattice and analyze the 
spatial and energetic distribution of their bound states. 
Details of the crystal growth, and the STM/STS methods 
have been described elsewhere^^. 

Fig. 1 a) shows an atomic resolution topographic im- 
age with two types of defects at the surface of LiFeAs. 
For a clear identification of the Fe, Li and As positions, 
and thereby assignment of the visible defect sites, a close 
look at the crystallography and the defects is needed. 

LiFeAs presents a layered crystal structure detailed in 
Fig. 1 . At its core is a square lattice of Fe atoms, 
each nested at the center of a tetrahedron of As. After 
cleaving between the weakly bonded Li layers, the surface 
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FIG. 1. Scanning tunneling topographies of typical defects 
found at the surface of LiFeAs. The estimated atomic con- 
figuration of the crystal structure with a — (3.74 zb 0.03) A^^ 
is shown on top of the topographies to locate the defect posi- 
tions (Li: yellow, As: blue, Fe: red), (a) 3 different Fe-D2-1 
and one Fe-C2 defect, (b) Schematic crystal structure in edge- 
on view of the upper layers of LiFeAs with Li at the top as 
expected from cleaving, (c), (d) and (e) are topographies of 
Fe-D2-2, As-Di and Li-Di defects. The tunneling conditions 
were It — 100 pA for all scans and Ub — 6 meV, z-range 120 
pm (a and e), Ub = — 12 meV, z-range 55 pm (c), Ub = — 12 
meV, z-range 120 pm (d). 



consists of a top layer of Li in a square lattice followed by 
a square lattice of As of the same periodicity but trans- 
lated by [1/2, 1/2]. The third layer is Fe arranged again 
in a square lattice but of twice the density and rotated 
45° relative to the Li or As lattice, where neither As nor 
Li sits directly above an Fe site. 

The most common defect (labeled Fe-Ds-l in Fig. [l^)), 
provides a key clue to the identification of the lattice 
observed by STM. Its distinctive feature is the dihedral 
D2 symmetry exhibiting two bright lobes oriented either 
along [100] or [010]. The symmetries discussed in this 
work are based on the SchonflieB notation adapted for 
two dimensions, as STM provides a weighted 2D pro- 
jection of the 3D crystal most sensitive to disruptions 
occurring on the upper rather than the buried planes. 

The question is: Which chemical site in the lattice can 
generate the observed topography? To form the center 
of this D2 defect, both the Li and As sites are unlikely 
as they present D4 symmetry due to the four nearest 
neighbors within each plane as discussed above. Iron 
sites, however, would induce a D2 -symmetry in the up- 
per planes due to the tetrahedrally bonded As nearest 
neighbors with two As in the plane above and two in the 
plane below (Fig. lb). All four As nearest neighbors 
would be affected by charge transfer, but STM will be 



most sensitive to changes on the upper two As. In the 
following we refer to this defect as Fe-D2-1 based on the 
apparent center position and symmetry. 

The spatial assignment of the Fe-D2-1 defect allows 
the registry of the three atomic sublattices (Fe, As and 
Li) to the measured STM corrugation. In Fig. la) the 
STM topography is overlaid with a model of the (001) 
cleaved crystal structure of LiFeAs (Li: yellow. As: blue, 
and Fe: orange), where the Fe site sits at the center of 
the Fe-D2-1 defect. This model works remarkably well 
across the full scan with no phase shift between the ap- 
plied grid and the measured atomic periodicity, matching 
the smaller Fe periodicity. Consequently, Fe-D2-1 always 
extends along the line of the two strongly polarizable^ 
As nearest neighbors lying above the Fe layer. 

Defects with similar topography and densities have 
been observed in STM studies of LiFeAs grown 
elsewhere^^ The simplest explanation for this com- 
monly observed defect are iron-site defects, even though 
other origins such as dimers or interstitial impurities can- 
not be completely excluded. Likely possibilities are Fe 
vacancies or Li substitution of iron sites, as the crystal 
was grown from a Li-rich flux. 

Previously, atomic resolution STM images of LiFeAs 
have produced the periodicity of either Li or A^^^H^. 
Here, the periodicity of the Fe lattice is observed at low 
bias voltages {\Ub\ ^ 20 meV) and with a particular tip 
wavefunction apparently yielding an enhanced sensitiv- 
ity to the iron in-plane d^^^-orbital. The latter becomes 
apparent in the measurement of the gaps by STS. A typ- 
ical STS spectrum taken at the center of a defect-free 
area is shown in Fig. 2a) in solid black. As reported 
previousl}^'^, two nodeless gaps are clearly resolved 
with half width peak-to-peak of A^^ = 6 meV and half 
width shoulder to shoulder = 3 meV, but the spec- 
tral weight contribution of the s mall g ap A2, originating 
from the in-plane Fe-d^j^y-orbitaJ^^^, is about twice as 
large as previously found^^. Further, the smaller atomic 
lattice corrugation of the Fe plane allows more accurate 
registry of topographic features with the crystal lattice. 

Turning to the effect of the Fe-D2-1 on the LDOS in 
the vicinity of the defect, the energetic and spatial dis- 
tribution of bound states are revealed. Fig. 2a) upper 
panel shows two raw spectra taken at 2.2 K at the center 
and on the lobes of the Fe-D2-1 defect as well as a ref- 
erence spectrum of pristine LiFeAs, taken from the same 
region. The lower panel presents these two spectra nor- 
malized by subtracting the reference to enhance subtle 
features^^. The first effect of the defect is a suppres- 
sion of spectral weight at the large gap coherence peaks. 
More interestingly, both spectra reveal a resonance at 3 
meV coinciding with the edge of the small gap Af^ and 
slightly more pronounced at the lobe than at the defect 
center. A 2 is on a Fermi surface segment dominated by 
the iron in-plane dxy band, expected to be influenced by 
defects in the iron plane. A bias symmetric counterpart 
at -3 meV is not visible within our resolution. The 4.2 K 
spatial distribution of the 3 meV bound state resonance. 



3 



(a)i5 

CO 10 
> 



Fe-D^-I 


1 / ' 











(b) Topography (a) 



10 




-10 



-5 



^2 

Sample Bias [meV] 



A^P 5 AfP 



10 



mm 



FIG. 2. a) The top panel shows dl/dV spectra taken at a 
Fe-D2-1 defect (red and blue lines) and the average over 20 
spectra measured about 2 nm away from the defect (black 
line). The thin red and blue lines in the lower panel corre- 
spond to the difference between the spectra taken at the de- 
fect and the average from the defect-free position. A bound 
state is visible at approximately 3 meV. Spectra were taken 
at the positions marked by single red and blue dots in the 
topographic image b). c) dl/dV map at an energy of 3 ± 1 
meV, i.e. averaged from 2 meV to 4 meV. The topography 
and the dl/dV maps are overlayed with the top- view crystal 
structure. Spectra shown in a) were taken at a temperature of 
2.2 K. Topography and dl/dV map shown in b) and c) were 
recorded at 4.2 K with a constant height defined by Ub,o = 25 
meV, Iq = 260 pA. Topography and dl/dV maps have edge 
lengths of ~ 34 nm. 



shown in Fig. 2 c), follows the D2 pattern measured in 
topography. To enhance spatial features, the maps were 
averaged over the 2 meV width of the bound state, i.e. 
from 2 meV to 4 meV for the 3 meV map. 

We now address four other recurrent defects of LiFeAs 
which we register based on our identification of the lat- 
tice as shown in Fig[l^), and c-e). A second Fe-centered 
defect also has a dihedral D2 symmetry and is referred to 
as Fe-D2-2 (Fig. lb). Fe-D2-2, which we will not discuss 
in detail, exhibits similar electronic properties as Fe-D2- 
1 but with a weaker 3 meV bound state resonance. The 
others are labeled Fe-C2, Li-Di and As-Di based on the 
apparent defect site registry in the xy-plane and their two 
dimensional group symmetry. We observed the two possi- 
ble chiralities and two orientations of Fe-C2 as well as all 
four possible orientations of Li-Di and As-Di. We note 
that except for the two Fe-D2 defects, the other three 
break the local symmetry of the crystal lattice. This 
symmetry breaking can be caused by particular orbital 
orientations close to the defect^^, by interstitial impurity 
positions or by the occurrence of two or more impurities 
close to each other. A summary of the defect character- 
istics is presented in Table-|l[ 

Fig 3a) shows spectra taken at the center and the arm 
of an As-Di defect. Similar to Fe-D2-1, As-Di causes a 
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FIG. 3. Same as Fig. [2] but for the As-Di defect. Bound 
states are visible at approximately -1.5 meV and at + 3.7 
meV. The topography which indicates the location of the 
spectra is shown as an inset of a), b) and c) dl/dV maps 
at energies of —1.5 and 3.7 ± 1 meV, respectively. 



strong suppression of the coherence peaks of Ai. How- 
ever, the spectra not only reveal a resonance close to the 
positive edge of A2 but also indicate a peak below -A2. 
Additionally, another resonance appears inside the gap 
at about -1.5 meV. The -1.5 meV resonance is present 
only in the arm of the defect but not at its center. This 
is made more clear by the spectral maps acquired at -1.5 
meV and 3.7 meV presented on Fig. 3b) and 3c). 
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FIG. 4. Same as Fig. [2]but for the Li-Di defect, a) A bound 
state is visible at approximately +1.2 meV. The topography 
which indicates the location of the spectra is shown as an 
inset of a), b), c) and d) show dl/dV maps at energies of -1, 
1 and 3.7 =b 1 meV, respectively. 

Spectra taken on the Li-Di defect are shown in Fig. 
4a). Again, the coherence peaks of the large gap are re- 
duced but a resonance near A2 is just weakly suggested. 
The three spectra in Fig. 4a) are qualitatively similar. 
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TABLE L Summary of the five defect types. The density is 
the number per LiFeAs formula unit. Eb is the bound state 
energy relative to Ep- 



Defect 


Symmetry 


Density 


Eb 


\Eb\ 


\Eb\ 






[lU J 


[meV] 


APP 
^1 


APP 


Fe-D2-1 


D2 


^ 1.2 ±0.3 


^ 3.0 


0.5 


1 


Fe-D2-2 


D2 


^ 0.2 ±0.1 


^ 3.0 


0.5 


1 


As-Di 


Di 


^ 0.5 ±0.1 


^ 3.7 


0.62 


1.23 








^-1.5 


0.25 


0.5 


Li-Di 


Di 


^ 0.1 ±0.1 


^ 3.7 


0.62 


1.23 








^ 1.2 


0.2 


0.4 


Fe-Cs 


C2 


^ 0.1 ±0.1 


^ ±3.5 


0.58 


1.17 








^ 1.0 


0.17 


0.33 



showing a clear in-gap resonance at about 1.2 meV. Spec- 
tral maps at 1 meV, at the bias symmetric energy -1 meV 
and at the 3.7 meV peak just above A2 are shown in Fig. 
4b), c) and d), revealing that the near 1 meV state is the 
most extended. At -1 meV, the LDOS is predominantly 
localized on two As sites near the Li center. 




FIG. 5. Same as Fig. [2] but for the Fe-C2 defect, a) Bound 
states are visible at approximately ±1.0 meV and at about ± 
3.5 meV. The topography which indicates the location of the 
spectra is shown as an inset of a), b), c) and d) show dl/dV 
maps at energies of -1, 1 and 3.7 ±1 meV, respectively. 

Fig. 5a) shows spectra recorded on three sites of a 
Fe-C2 defect. Once more, the defect reduces the ampli- 
tude of the coherence peaks of Ai. Resonances at 1 meV 
and approximately ±3.5 meV are observed. The spa- 
tial dl/dV mapping of these defects differs from the con- 
stant current map topography. The dl/dV map at the 
-3.5 meV resonance (not shown) is qualitatively similar 
to the -1 meV map. Interestingly, the —1 meV and 3.7 
meV maps show one chirality while the opposite chirality 
is observed at ±1 meV. 

Summarizing the properties of all defects, some strik- 
ing features of the dl/dV maps presented here should be 
highlighted. Bound state resonances exhibit large spa- 



tial extents, most distinct within an area of ~ 4 x 4 
unit cells. Three of the five defect types observed in 
LiFeAs break the lattice symmetry and the bound state 
patterns are not commensurate with the crystal lattice. 
The strong dependence on the bias voltage reveals dif- 
ferent spatial localization of different resonance energies. 
Consequently, bound states in LiFeAs result in more com- 
plex patterns than those observed in the cuprates where 
the LDOS of particle and hole components shows com- 
plementary patterns overlaying the CuO lattice site^. 
In particular, the chiral C2 pattern of the Fe-C2 defect 
is rarely observed. A recent observation of this in NbSe2 
is attributed to charge density waves^^, and it has been 
recently proposed that impurities in pnictides may yield 
chiral C2 LDOS patterns in the presence of spin density 
orde^^ or strong orbital fluctuations^^. If the extraor- 
dinary chiral pattern of the Fe-C2 defect is indeed con- 
nected to magnetic, charge or orbital order locally frozen 
by the defect, the present observations might provide fur- 
ther insight regarding the pairing glue in LiFeAs. 

We turn to the analysis of bound state energies rela- 
tive to the gap sizes, which can help to identify the gap 
structure. Both of the two lattice symmetry preserving 
Fe-D2 defects stand out as having only one bound state 
appearing near the edge of the small gap. AS-D2, Li- 
Di and Fe-C2 also cause resonances at or near the small 
gap edge A2, but lower temperature measurements are 
needed to clearly resolve these features. AS-D2, Li-Di 
and Fe-C2, which break the local lattice symmetry, in- 
duce pair-breaking indicated by in-gap bound states pro- 
nounced at either positive or negative bias. For most 
resonances only one peak is visible and the expected 
bias symmetric counterpart^ ^ is either not observed or 
only weakly present, revealing a significant particle- hole 
asymmetry of the bound states as observed by STS. 
Due to our resolution and the uncertainty inherent in 
finding the bound state peaks on a background that is 
also influenced by the presence of the defect, we cannot 
definitively ascribe these states to the small gap edge, 
and cannot exclude that they are in-gap states of the 
larger gap. However, the asymmetric multi-peak ener- 
getic structure we observe and the observation of states 
near the small gap are consistent with recent theoretical 
predictions for s^ symmetries'^. While we cannot ex- 
clude the possibility that the defects in the crystal exhibit 
magnetic properties, potentially inducing in-gap bound 
states for either s^ and s+^l^, such states at the gap 
edge have not previous ly been observed in conventional 
s-wave super conduct or ^^SElM]^ Controlled introduction 
of impurities and probing of the bound state resonances 
should provide further evidence for the pairing symmetry 
in the iron pnictides. 

Regardless of pairing symmetry and interaction, we 
gain information regarding the resilience of the super- 
conductive phase of LiFeAs to chemical substitution. In 
analogy to the observation of pair-breaking and non-pair 
breaking defects in BaFe2As2, ^^^' *the two Fe-D2 defects 
seem to correspond to such non-pair breaking impurities. 
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at least for the small gap, while the Fe-C2 defect forms 
a pair-breaking substituent also likely located within the 
iron layer. 

In conclusion, we classify and register the common de- 
fects appearing at the surface of LiFeAs. Our registry 
can be used as a guideline for the identification of delib- 
erately introduced external impurities. The bound states 
show complex spatial patterns that vary strongly with en- 
ergy and do not reflect the lattice periodicity. Further, 
the energetic multi-peak bound state structures with res- 
onances at the small superconducting gap are in accor- 



dance with calculations of non-magnetic defects in s^ gap 
structures^, although the possibility of a magnetic nature 
of our defects does not allow a definitive attribution of 
the gap structure. Finally, the observed variety of dif- 
ferent defects suggests many possibilities for controlled 
doping and thereby for tailoring the material properties. 
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